• 专利附录
  • 专利说明
  • 权利要求
  • 类似技术
  • 同族专利
  • 引用文献
  • 法律状态
  • 优先权
    • 专利摘要:
    La présente invention concerne une méthode de détection pour un récepteur de système MIMO dans lequel on effectue une détection linéaire (ZF,MMSE,MF) pour fournir un vecteur égalisé. Ce vecteur égalisé est représenté dans une base réduite obtenue à partir de la réduction de la matrice de canal. Il fait l'objet d'un processus itératif de suppression de bruit dans la représentation selon la base réduite. A chaque itération, on recherche la composante du vecteur égalisé dans la base réduite la plus éloignée d'une zone non bruitée (Z<sub>ε</sub>) entourant la constellation produit avec une marge de tolérance, et le point représentatif du vecteur égalisé de cette zone en lui soustrayant un vecteur de bruit dans la direction de cette composante, dont le module est égal à une fraction de la marge de tolérance. La suppression itérative converge lorsque le vecteur égalisé appartient à la zone non bruitée. Le vecteur égalisé ainsi obtenu peut faire l'objet d'une décision dure par quantification sur les points de la constellation produit ou bien d'une décision douce pour calculer les valeurs souples des bits ayant servi à générer le symbole transmis.
    公开号:EP3696992A1
    申请号:EP20151153.2
    申请日:2020-01-10
    公开日:2020-08-19
    发明作者:David Demmer；Jean-Baptiste Dore
    申请人:Commissariat a lEnergie Atomique et aux Energies Alternatives；
    IPC主号:H04B7-00
    专利说明:
    [0001] La présente invention concerne le domaine des systèmes de télécommunication multi-antenne encore dénommés MIMO (Multiple Input Multiple Output) et plus particulièrement ceux utilisant une détection aidée par réduction de réseau (Lattice Reduction aided MIMO detection). ÉTAT DE LA TECHNIQUE ANTÉRIEURE
    [0002] Les systèmes de télécommunication multi-antenne ou MIMO (Multiple Input Multiple Output) sont bien connus de l'état de la technique et ont fait leur entrée, il y a une dizaine d'années, dans plusieurs normes de systèmes de télécommunication sans fil (par exemple IEEE 802.11n). Ces systèmes présentent, en raison de leur capacité de multiplexage et leur diversité spatiale, des capacités de canal sensiblement supérieures aux systèmes conventionnels mono-antenne SISO (Single Input Single Output).
    [0003] Un système de transmission MIMO à NT antennes d'émission et NR antennes de réception est caractérisé par sa matrice de canal H, de taille NR ×NT dont les éléments sont les coefficients d'atténuation complexes des canaux élémentaires entre les NT antennes d'émission et NR antennes de réception.
    [0004] On a représenté de manière schématique en Fig. 1 un tel système MIMO. Un élément hik de la matrice H donne le coefficient de canal élémentaire entre l'antenne d'émission d'indice k et l'antenne de réception d'indice i.
    [0005] La transmission d'un bloc de symboles dans un tel système peut être modélisée par :y = Hs + n
    [0006] Les symboles émis appartiennent à un alphabet de modulation, noté A, par exemple un alphabet M -QAM.
    [0007] On supposera dans la suite que le nombre d'antennes de réception est supérieur ou égal au nombre d'antennes d'émission, c'est-à-dire NR ≥ NT.
    [0008] Les éléments des vecteurs et de la matrice intervenant dans la relation (1) étant des valeurs complexes, on traite leurs parties réelles et imaginaires séparément. Autrement dit, si y , s et n sont les vecteurs réels respectivement de taille nR, nT et nR avec nR = 2NR et nT = 2NT définis par :
    [0009] Différents types de détection peuvent être prévus à la réception.
    [0010] La détection optimale utilise un critère de maximum de vraisemblance ou ML (Maximum Likelihood), consistant à rechercher sur l'ensemble des vecteurs de symboles émis possibles celui minimisant la distance euclidienne au vecteur reçu, autrement dit : s ^ML = arg min s ∈ A N T‖ y − Hs ‖
    [0011] Toutefois, cette recherche présente une complexité en MnT et n'est donc pas réalisable en pratique pour des ordres de modulation élevés et un grand nombre d'antennes d'émission.
    [0012] En pratique, on a recours à des méthodes de détection sous-optimales, dites linéaires, permettant de séparer les différents symboles au moyen de combinaisons linéaires. Plus précisément, revenant à la formulation complexe de la relation (1), le vecteur des symboles émis est estimé en multipliant le vecteur des signaux reçus par une matrice de détection (encore dénommée matrice d'égalisation) W de taille nT × nR :r = Wy = WHs + Wn
    [0013] Les avantages et inconvénients respectifs de ces différents types de détection sont connus : ainsi une détection de type ZF peut amplifier considérablement les composantes de bruit et une détection de type MMSE, si elle réduit cette amplification, nécessite l'estimation du rapport signal sur bruit et est donc plus complexe.
    [0014] Les symboles émis sont ensuite estimés en prenant une décision (dec) dure ou douce sur les composantes du vecteur égalisé, r : s ^ZF = dec W ZF y
    [0015] Par exemple dans le cas d'une décision dure, pour chaque composante du vecteur égalisé, on recherche l'élément de l'alphabet, c'est-à-dire le point de la constellation de modulation le plus proche de cette composante.
    [0016] Si un détecteur linéaire ou LD (Linear Detector) est relativement simple, ses performances sont fonction du conditionnement de la matrice H et plus précisément son écart à l'orthogonalité. En effet, si les valeurs propres de la matrice H H H sont faibles les composantes de bruit (Wn dans l'expression (3)) sont amplifiées et le taux d'erreur (BER) est élevé.
    [0017] Afin de remédier à cette difficulté, il a été proposé d'effectuer une détection linéaire aidée par réduction de réseau de points comme décrit dans l'article de D. Wübben intitulé « Near-maximum likelihood détection of MIMO systems using MMSE-based lattice reduction » publié dans Proc. IEEE International Conférence on Communications, vol. 2, Juin 2004, pp. 798-802.
    [0018] On rappelle tout d'abord qu'un réseau de points A dans
    [0019] La base V génératrice du réseau est représentée par la matrice V = (v 1|v 2|...|v n ) dans la base canonique de
    [0020] Un même réseau de points, Λ, peut être généré par des bases différentes. Ainsi, si V est une base génératrice du réseau, représentée par la matrice V et si T est une matrice unimodulaire (autrement dit une matrice à coefficients entiers et dont la valeur absolue du déterminant est égale à 1), les colonnes de la matrice V' définie par :Vʹ = VT
    [0021] La réduction d'un réseau de points Λ consiste à trouver une base génératrice de ce réseau formée par des vecteurs de plus faible longueur et les plus proches de l'orthogonalité.
    [0022] L'écart à l'orthogonalité de la base V est donné par l'indicateur :od V = ∏ i = 1 n ‖ v i ‖ vol Λ
    [0023] La réduction du réseau de points peut être approchée au moyen de l'algorithme LLL(Lenstra-Lenstra-Lovász). Une base du réseau LLL-réduite, Ṽ de paramètre δ du réseau est obtenue à partir d'une décomposition QR de la matrice V . Cette base est donnée par les vecteurs colonnes de la matrice Ṽ = Q̃R̃ vérifiant :r ˜ l , k≤ 1 2 r ˜ l , l pour 1 ≤ l ≤ k ≤ n
    [0024] Le paramètre δ donne la qualité de la base Ṽ, l'approximation est d'autant meilleure que la valeur δ est proche de 1.
    [0025] Les points d'une constellation de modulation M -QAM, peuvent être considérés, à une translation et un facteur multiplicatif communs près, comme des éléments de
    [0026] La réduction du réseau et notamment sa réduction LLL permet d'obtenir une base de ce réseau composée de vecteurs courts et proches de la condition d'orthogonalité. Si H̃ est la matrice du canal représenté dans cette nouvelle base, on a : H ˜ = H ‾ T
    [0027] La relation (4) peut alors se ré-écrire sous la forme : y ‾ = H ‾ T T − 1s ‾+ n ‾ = H ˜z ‾ + n ‾
    [0028] Si l'on effectue une détection linéaire à la réception, on a dans le domaine réel : r ‾ = W ‾y ‾ = W ‾H ˜z ‾ + W ‾n ‾
    [0029] Pour des raisons de simplification des notations, nous abandonnerons dans la suite les traits suscrits pour les vecteurs et matrices réelles. Ainsi, la relation (15) s'écrit sous forme simplifiée :r = W ˜ y = W ˜H ˜ z + W ˜ n
    [0030] L'homme du métier comprendra toutefois que les vecteurs et matrices apparaissant dans l'équation (16) sont à valeurs réelles et que leurs dimensions sont doubles de celles figurant dans l'équation (6).
    [0031] Une première méthode de détection linéaire aidée par réduction de réseau de points ou LR-LD (Lattice Reduction aided Linear Detection) a été décrite dans l'article précité. Elle consiste à quantifier chaque composante zi, i = 1,...,nT de z sur le ième vecteur de la base réduite. La quantification de la composante zi est réalisée en recherchant le point appartenant à la projection de la constellation produit ANT sur le i-ème vecteur. En d'autres termes, dans chacune des directions portées par un vecteur de la base réduite, on recherche le point de la projection de la constellation ANT sur cette direction qui est le plus proche de la composante de z sur ce vecteur.
    [0032] Une seconde méthode de détection linéaire aidée par réduction de réseau, également décrite dans l'article précité, dénommée ci-après LR-LD-SIC (Lattice Reduced Linear Détection with Successive Interference Cancellation) consiste à effectuer une décomposition QR de la matrice de canal H̃ (où Q ,R sont respectivement une matrice orthogonale et une matrice triangulaire supérieure) puis de détecter les composantes de z, en commençant par la composante znT non interférée par les autres composantes (du fait de la structure triangulaire supérieure de R). Les colonnes de H̃ peuvent être permutées (ce qui est équivalent à une permutation des antennes d'émission) de manière à ce que le rapport signal sur bruit de la composante znT soit le plus élevé possible. Autrement dit, on effectue pour chaque permutation des colonnes de H̃ une décomposition QR et l'on sélectionne la permutation conduisant au coefficient rnTnT le plus élevé.
    [0033] La méthode LR-LD présente l'inconvénient d'opérer une troncation sous-optimale dans la mesure où selon l'ordre selon lequel les composantes de z sont tronquées, les performances en termes de BER peuvent être très différentes, ce pour une même valeur de rapport signal sur bruit (SNR).
    [0034] La méthode LR-LD-SIC est performante lorsque la détection linéaire est de type MMSE mais elle nécessite une importante quantité de calculs, en particulier pour des configurations MIMO ayant un grand nombre d'antennes, ce qui n'est guère compatible avec la capacité calculatoire d'un terminal mobile.
    [0035] En outre, ni la méthode LR-LD ni la méthode LR-LD-SD ne sont à même de fournir des sorties souples, telles que celles nécessaires à décodage canal à entrées souples lorsqu'un turbocode, un code convolutif ou un code LDPC est utilisé à l'émission.
    [0036] Un but de la présente invention est par conséquent de proposer une méthode de détection pour un récepteur de système MIMO aidée par réduction de réseau de points, qui présente de bonnes performances en termes de taux d'erreur en fonction de rapport signal sur bruit, sans nécessiter de calculs importants. Un but subsidiaire de la présente invention est de proposer une telle méthode de détection MIMO aidée par réduction de réseau de points qui soit en outre capable de fournir un résultat de détection sous forme de valeurs souples. EXPOSÉ DE L'INVENTION
    [0037] La présente invention est définie par une méthode de détection pour récepteur de système MIMO à NT antennes d'émission et NR antennes de réception, caractérisé par une matrice de canal H, chaque antenne d'émission transmettant un symbole appartenant à une constellation de modulation, le vecteur des symboles émis appartenant à une constellation produit dans l'espace
    [0038] La matrice de détection linéaire est notamment choisie parmi des matrices de détection ZF, MMSE ou MF définies respectivement par : W ˜ZF =H ˜† = H ˜H H ˜ − 1 H ˜H
    [0039] Les constellations de modulation relatives aux différentes antennes d'émission sont typiquement des constellations QAM. Les constellations de modulation relatives aux différentes antennes sont avantageusement choisies identiques.
    [0040] Avantageusement, pour déterminer la composante du vecteur égalisé dans la base réduite, la plus éloignée de la zone non bruitée, on recherche l'indice ℓ tel queℓ = arg max ic i − s max − ε / ν i
    [0041] Le vecteur de bruit soustrait au vecteur égalisé c peut être déterminé par a ℓ sgn c ℓ. μ . ε v ℓ
    [0042] Selon une première variante, on effectue une décision dure sur le vecteur égalisé au terme du processus itératif, en recherchant le point de la constellation produit le plus proche de ce vecteur.
    [0043] Selon une seconde variante, préalablement au processus itératif on initialise les variances, σ n , m 2 ,
    [0044] Dans ce dernier cas, on peut calculer pour chaque symbole émis par une antenne d'émission m les valeurs souples des différents bits bm,q, bm,q étant le qème bit du mot ayant servi à générer le symbole émis s m , les valeurs souples étant données parL b m , q= 1 2 σ n , m 2 min u ∈ S m , q 1 ‖ c − u ‖ 2 − min u ∈ S m , q 0 ‖ c − u ‖ 2
    [0045] D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture d'un mode de réalisation préférentiel de l'invention, décrit en référence aux figures jointes parmi lesquelles : La Fig. 1 représente de manière schématique un système MIMO connu de l'état de la technique ; La Fig. 2 représente de manière schématique un ordinogramme d'une méthode de détection linéaire aidé par réseau de points et réduction successive de bruit selon un premier mode de réalisation de l'invention ; La Fig. 3 illustre de manière schématique un exemple de réduction de bruit dans une étape la méthode de la Fig. 2 ; La Fig. 4 représente de manière schématique un ordinogramme d'une méthode de détection linéaire aidé par réseau de points et réduction successive de bruit selon un second mode de réalisation de l'invention ; La Fig. 5 représente le débit du système MIMO en fonction du rapport signal sur bruit pour un procédé de détection selon le second mode de réalisation de l'invention et selon un procédé de détection selon l'état de la technique. EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS
    [0046] On considérera dans la suite un système MIMO comme décrit en relation avec la Fig. 1. Les symboles d'information émis par les NT antennes d'émission appartiennent à des constellations QAM. Sans perte de généralité, on supposera que l'alphabet (ou constellation) de modulation utilisé(e) pour les différentes antennes d'émission est identique, en d'autres termes un bloc de symboles émis appartient à une constellation produit ANT . On supposera également que la puissance d'émission de chaque antenne est normalisée à 1. En absence de bruit, le signal reçu peut être représenté par un point du réseau de points Λ H engendré par les vecteurs-colonnes de H dans l'espace
    [0047] On supposera en outre, comme précédemment, que les éléments de la matrice H sont connus du récepteur, par exemple au moyen d'une estimation de canal réalisée à partir de symboles pilotes.
    [0048] L'idée à la base de l'invention est de procéder à une suppression de bruit par réduction itérative du bruit sur les différentes composantes du signal égalisé, dans la base réduite du réseau de points.
    [0049] La Fig. 2 représente de manière schématique un ordinogramme d'une méthode de détection linéaire aidée par réseau de points et réduction successive de bruit selon un premier mode de réalisation de l'invention.
    [0050] A l'étape 210, on effectue une réduction LLL du réseau de points Λ H engendré par les vecteurs colonnes de la matrice de canal, exprimée sous forme réelle, c'est-à-dire de taille nR ×nT . Cette réduction LLL est connue en soi, elle utilise une décomposition QR de la matrice de canal (H = QR où Q est une matrice orthogonale et R une matrice triangulaire supérieure) et fournit en retour des matrices Q̃, R̃ et T avec H̃ = Q̃R̃ où H̃ est la base (δ - )LLL réduite de Λ H (ou matrice de canal dans la base réduite), Q̃ est une matrice orthogonale et R̃ une matrice triangulaire supérieure et T est une matrice de transformation unimodulaire de taille nT × nT telle que H̃ = HT, les vecteurs colonnes de T formant la base réduite.
    [0051] A l'étape 220, on effectue une détection linéaire du signal reçu au moyen d'une matrice d'égalisation W̃.
    [0052] Cette détection linéaire pourra être par exemple une détection de type ZF ( W̃ = W̃ ZF ), de type MMSE (W̃ = W̃ MMSE ) ou de type MF (W̃ = W̃ MF ), de manière connue en soi. Si l'on note y le vecteur reçu par les antennes de réception et r = W̃y le vecteur égalisé. Ce dernier peut être considéré comme un vecteur réel à nT = 2NT éléments.
    [0053] A l'étape 230, on obtient une représentation du signal égalisé dans la base réduite en calculant le vecteur c = Tr de taille nT .
    [0054] A l'étape 240, on calcule la matrice de Gram A = TT T de la base réduite et si la matrice de Gram est non-diagonale, on passe à l'étape 250. La matrice de Gram donne les produits scalaires des vecteurs colonnes de T .
    [0055] On notera que si la matrice A est diagonale, les vecteurs colonnes de T sont orthogonaux entre eux. Dans ce cas, on termine la détection à l'étape 290.
    [0056] A l'étape 250, on détermine pour chaque vecteur de la base réduite, c'est-à-dire chaque vecteur colonne t i , i = 1,...,nT, de la matrice T le vecteur colonne t k tel que :k = arg max j ≠ i 〈 t i | t j 〉
    [0057] On entre ensuite dans une boucle itérative de réduction de bruit en 260.
    [0058] A l'étape 260, on détermine si la condition suivante est vérifiée :∃ i ∈ 1 , … , n T tel que c i > s max + ε / ν i
    [0059] Si la condition (18) n'est pas vérifiée, autrement dit si toutes les composantes du vecteur égalisé vérifient ci ≤ s max + ε/vi on termine la détection en 290.
    [0060] A défaut, on détermine en 270 la composante ℓ du vecteur égalisé qui s'écarte le plus de la constellation produit, plus précisément qui s'écarte le plus d'une zone de sécurité (dite aussi de réception non bruitée) Zε dans l'espace
    [0061] En d'autres termes, on recherche la direction d'erreur maximale par rapport à la constellation produit. Il est essentiel de noter que cette recherche est effectuée dans la base réduite, autrement dit, dans une base non-orthogonale.
    [0062] On effectue ensuite en 280 une réduction de bruit en soustrayant au vecteur égalisé une fraction µ de la marge d'erreur selon la direction de la composante déterminée à l'étape précédente :c = c − a ℓ sgn c ℓ. μ ε v ℓ
    [0063] Autrement dit, on rapproche le vecteur égalisé de la zone non bruitée Z ε en réduisant la distance à cette zone dans la direction d'erreur maximale. On utilise de préférence une fraction de 1 2
    [0064] On notera que le vecteur égalisé, mis à jour au moyen de la relation (20), ne voit pas seulement sa composante c ℓ être modifiée puisque A n'est pas une matrice diagonale.
    [0065] Ainsi l'algorithme de détection pourra passer plusieurs fois sur la même direction, en corrigeant une erreur introduite sur une autre direction.
    [0066] L'algorithme se poursuit en retournant à l'étape de test 260.
    [0067] Les itérations se terminent lorsque le vecteur égalisé se trouve dans la zone de réception non bruitée Zε (auquel cas le bruit aura été supprimé) ou, à défaut, lorsqu'un nombre maximal d'itérations est atteint (absence de convergence). En cas de convergence, le vecteur égalisé se trouve au pire à une distance (1 - µ)ε de la constellation et dans le cas oùμ = 1 2 ,
    [0068] A l'étape 290, on effectue une quantification du vecteur égalisé ayant fait l'objet de la réduction itérative de bruit. Cette quantification est réalisée comme dans l'art antérieur, en recherchant le point de la constellation le plus proche du point représentant le vecteur égalisé.
    [0069] La Fig. 3 illustre sur un exemple l'étape de réduction de bruit 280. Dans cet exemple, on a supposé que l'on avait une seule antenne d'émission (nT = 2) et que la constellation de modulation utilisée était une constellation 4-QAM. La constellation produit se réduit donc à cette constellation 4-QAM.
    [0070] La constellation produit a été représentée dans la base réduite, autrement dit la figure représente l'image de l'alphabet de modulation A par la transformation T -1.
    [0071] La base réduite est constituée des vecteurs t 1 = 0 1
    [0072] Les symboles de modulation non bruités, c'est-à-dire les points de la constellation de modulation sont représentés en A, B, C, D. Le signal reçu, après égalisation, c'est-à-dire le vecteur r, est représenté par le point M. Le vecteur c est supposé être ici− 2 − 1 .
    [0073] La zone de réception non bruitée, autrement dit la zone Zε entoure la constellation de modulation avec une marge d'erreur ε.
    [0074] Selon la direction du vecteur t 2, le point M appartient à la zone Zε et donc on ne réduit pas le bruit selon cette direction.
    [0075] En revanche, selon la direction du vecteur t 1, le point M n'appartient pas à à la zone Zε et donc il est ramené vers cette zone selon cette direction au moyen de l'expression (16) oùμ = 1 2 .
    [0076] On voit donc que le nouveau vecteur égalisé est plus proche de la zone Zε de ε 2
    [0077] Dans un second mode de réalisation de l'invention, on procède à une détection souple des symboles émis. Ce sera par exemple le cas lorsque les symboles successifs émis sont le résultat d'un codage d'une séquence binaire, par exemple du turbocodage d'un paquet de bits. Dans un tel cas, le décodeur convolutif (utilisant par exemple un algorithme SOVA ou un algorithme de Viterbi classique) nécessite de recevoir en entrée les valeurs souples des bits de la séquence compte tenu des symboles reçus.
    [0078] On se place tout d'abord dans le cas d'un système SISO, c'est-à-dire à une seule antenne d'émission et une seule antenne de réception.
    [0079] Si l'on note bm,q le qème bit du mot ayant servi à générer le symbole émis s m appartenant à l'alphabet de modulation A, le récepteur doit évaluer la valeur souple ou LLR (Log Likelihood Ratio), L(bm,q ) définie par :L b m , q= ln P b m , q = 1 | yP b m , q = 0 | y= ln P y | b m , q = 1P y | b m , q = 0
    [0080] Les probabilités conditionnelles P( y |bm,q = X) (avec X = 0 ou 1) sont déterminées à partir de :P y | b m , q = X = ∑ u ∈ S m , q X P y | s = u
    [0081] De manière classique, on peut approximer la somme des termes apparaissant à droite de la relation (23) par le terme prédominant (approximation dite Max-Log) :P y | b m , q = X ≈ max u ∈ S m , q X P y | s = u
    [0082] Compte tenu de la densité gaussienne de la probabilité conditionnelle, la valeur souple du bit bm,q s'exprime sous la forme suivante :L b m , q= 1 2 σ n 2 min u ∈ S m , q 1 ‖ y − Hu ‖ 2 − min u ∈ S m , q 0 ‖ y − Hu ‖ 2
    [0083] La relation (25) s'étend sans peine à un système MIMO utilisant une détection linéaire :L b m , q= 1 2 σ n , m 2 min u ∈ S m , q 1 ‖ r − u ‖ 2 − min u ∈ S m , q 0 ‖ r − u ‖ 2
    [0084] Le calcul des valeurs souples selon (26) et (27) doit être adapté dans le cas où on effectue une réduction itérative de bruit en sus de l'égalisation, comme illustré en relation avec la Fig. 4.
    [0085] La Fig. 4 représente de manière schématique un ordinogramme d'une méthode de détection linéaire aidé par réseau de points et réduction successive de bruit selon un second mode de réalisation de l'invention.
    [0086] Les étapes 410 à 480 sont identiques aux étapes 210 à 280 déjà décrites en relation avec la Fig. 2 et leur description ne sera donc pas reprise ici.
    [0087] En revanche, une étape d'initialisation a été ajoutée en 455 ainsi qu'une étape de mise à jour de la puissance de bruit en 475.
    [0088] Plus précisément, on initialise en 455 les variances de bruit (ou puissances de bruit) affectant les différentes composantes du vecteur égalisé au moyen de la relation (27).
    [0089] Ensuite, à chaque itération on met à jour en 475 les variances de bruit σ n , m 2
    [0090] Enfin à l'étape 490, au lieu de prendre une décision dure comme à l'étape 290, on calcule les distances du vecteur égalisé c aux différents points de la constellation produit puis l'on calcule des valeurs souples des différents bits au moyen de :L b m , q= 1 2 σ n , m 2 min u ∈ S m , q 1 ‖ c − u ‖ 2 − min u ∈ S m , q 0 ‖ c − u ‖ 2
    [0091] Si un turbocodage a été réalisé à l'émission pour générer des symboles sur chacune des antennes, on pourra procéder à un décodage de type BCJR utilisant en entrée les valeurs souples fournies par l'étape 490.
    [0092] La Fig. 5 représente le débit du système MIMO en fonction du rapport signal sur bruit pour un procédé de détection selon le second mode de réalisation de l'invention et selon un procédé de détection selon l'état de la technique.
    [0093] Le système MIMO utilisé était un système 2x2 avec une modulation 16-QAM. Un turbocodage de rendementR = 1 2
    [0094] Le débit du système est donné par ρ(1-η) Q où ρ est la vitesse de transmission en bit/s, Q est la taille du paquet et η est le taux d'erreur binaire.
    [0095] La courbe 510 représente le débit du système lorsque l'on effectue à la réception une détection linéaire de type ZF et un calcul de valeurs souples selon (22) et (23), c'est-à-dire sans l'approximation Max-Log, avant d'effectuer un décodage BCJR.
    [0096] La courbe 520 représente le débit du système lorsque l'on effectue à la réception une détection linéaire de type ZF, une suppression itérative de bruit selon le second mode de réalisation de l'invention, et un calcul des valeurs souples selon (22) et (23) avant d'effectuer un décodage BCJR.
    [0097] La courbe 530 représente le débit du système dans les mêmes conditions que 510 mais en utilisant l'approximation Max-Log, autrement dit en utilisant les expressions de valeurs souples selon (26) et (27).
    [0098] La courbe 540 représente le débit du système selon le second mode de réalisation de l'invention, mais en utilisant l'approximation Max-Log, autrement dit en utilisant les expressions de valeurs souples selon (26) et (27) et la mise à jour itérative des variances selon (28), avant d'effectuer un décodage canal.
    [0099] On remarque une amélioration importante du débit dans le cas où on utilise le second mode de réalisation de l'invention, cf. courbes (520) et (540), au lieu d'une méthode de détection de l'état de la technique, cf. courbes (510) et (530).
    权利要求:
    Claims (10)
    [0001] Méthode de détection pour récepteur de système MIMO à NT antennes d'émission et NR antennes de réception, caractérisé par une matrice de canal H, chaque antenne d'émission transmettant un symbole appartenant à une constellation de modulation, le vecteur des symboles émis appartenant à une constellation produit dans l'espace
    [0002] Méthode de détection selon la revendication 1, caractérisée en ce que la matrice de détection linéaire est choisie parmi des matrices de détection ZF, MMSE ou MF définies respectivement par : W ˜ZF =H ˜† = H ˜H H ˜ − 1 H ˜H
    [0003] Méthode de détection selon la revendication 1 ou 2, caractérisée en ce que les constellations de modulation relatives aux différentes antennes d'émission sont des constellations QAM.
    [0004] Méthode de détection selon la revendication 3, caractérisée en ce que les constellations de modulation relatives aux différentes antennes sont identiques.
    [0005] Méthode de détection selon la revendication 4, caractérisée en ce que pour déterminer la composante du vecteur égalisé dans la base réduite, la plus éloignée de la zone non bruitée, on recherche l'indice ℓ tel queℓ = arg max ic i − s max − ε / ν i
    [0006] Méthode de détection selon la revendication 5, caractérisée en ce que le vecteur de bruit soustrait au vecteur égalisé c est déterminé par a ℓ sgn c ℓ. μ . ε v ℓ
    [0007] Méthode de détection selon la revendication 6, caractérisée en ce queμ = 1 2 .
    [0008] Méthode de détection selon l'une des revendications précédentes, caractérisée en ce que l'on effectue une décision dure sur le vecteur égalisé au terme du processus itératif, en recherchant le point de la constellation produit le plus proche de ce vecteur.
    [0009] Méthode de détection selon la revendication 6, caractérisée en ce que préalablement au processus itératif, on initialise les variances, σ n , m 2 ,
    [0010] Méthode de détection selon la revendication 9, caractérisée en ce que l'on calcule pour chaque symbole émis par une antenne d'émission m les valeurs souples des différents bits bm,q, bm,q étant le qème bit du mot ayant servi à générer le symbole émis s m , les valeurs souples étant données parL b m , q= 1 2 σ n , m 2 min u ∈ S m , q 1 ‖ c − u ‖ 2 − min u ∈ S m , q 0 ‖ c − u ‖ 2 où σ n , m 2 est la valeur de la variance affectant la m-ème composante du vecteur égalisé c au terme du processus itératif, S m , q 0 est l'ensemble des points de la constellation produit admettant bm,q =0 et S m , q 1 est l'ensemble des points de la constellation produit admettant bm,q =1.
    类似技术:
    公开号 | 公开日 | 专利标题
    Wei et al.2016|Low complexity iterative receiver design for sparse code multiple access
    US9225451B2|2015-12-29|Transmission apparatus, reception apparatus, transmission method, reception method, and method for generating multi-dimensional constellations
    Steingrimsson et al.2003|Soft quasi-maximum-likelihood detection for multiple-antenna wireless channels
    US8711989B1|2014-04-29|Method and apparatus for receiving signals in a MIMO system with multiple channel encoders
    Ylioinas et al.2008|Iterative joint detection, decoding, and channel estimation in turbo-coded MIMO-OFDM
    JP6410948B2|2018-10-24|光信号を復号化する方法および受信機
    RU2414062C2|2011-03-10|Детектирование и декодирование с уменьшенной сложностью для приемника в системе связи
    EP2795826B1|2020-04-22|Procédé de transmission d&#39;un signal numérique pour un système ms-marc semi-orthogonal, produit programme et dispositif relais correspondants
    EP1547289B1|2009-04-29|Reception d&#39;un signal module selon une technique de codage multi-niveaux
    US7616699B2|2009-11-10|Method of soft bit metric calculation with direct matrix inversion MIMO detection
    EP1233565B1|2006-08-23|Procédé et appareil de réception de turbo-code pour un système MIMO
    US8279954B2|2012-10-02|Adaptive forward-backward soft output M-algorithm receiver structures
    US8867649B2|2014-10-21|Antenna layer selection and soft demapping for MIMO decoder
    JP5436453B2|2014-03-05|ツリーポジション適応の軟出力ｍアルゴリズムの受信機構造
    EP1212872B1|2008-12-17|Recepteur a multiplexage par repartition en frequences orthogonales avec estimation iterative de canal et procede correspondant
    EP1677434B1|2017-05-10|Détecteur à complexité réduite pour des systèmes à antennes multiples
    US8432987B2|2013-04-30|Method and apparatus for multiple antenna communications, and related systems and computer program
    US8000422B2|2011-08-16|Apparatus and method for detecting signal in multiple-input multiple-output | wireless communication system
    EP0654910B1|1998-08-19|Procédé de décodage itératif de codes en blocs concaténés
    KR100866805B1|2008-11-04|다중 안테나 통신시스템에서 로그 우도 율 생성 장치 및방법
    KR101124863B1|2012-03-27|다중 소스로부터의 통신신호를 처리하는 장치 및 방법
    EP1495608B1|2012-08-08|Architecture de decodeur de type egaliseur/treillis a concatenation pour un recepteur de television numerique
    US8184993B2|2012-05-22|Polarization mode dispersion | compensation in polarization multiplexed coded orthogonal frequency division multiplexing | systems
    Baek et al.2009|Combined QRD-M and DFE detection technique for simple and efficient signal detection in MIMO-OFDM systems
    EP1691492A1|2006-08-16|Sélection du détecteur en fonction de la qualité du signal reçu
    同族专利:
    公开号 | 公开日
    FR3092716A1|2020-08-14|
    FR3092716B1|2021-01-22|
    US20200259544A1|2020-08-13|
    引用文献:
    公开号 | 申请日 | 公开日 | 申请人 | 专利标题
    法律状态:
    2020-07-17| PUAI| Public reference made under article 153(3) epc to a published international application that has entered the european phase|Free format text: ORIGINAL CODE: 0009012 |
    2020-07-17| STAA| Information on the status of an ep patent application or granted ep patent|Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE |
    2020-08-19| 17P| Request for examination filed|Effective date: 20200110 |
    2020-08-19| AK| Designated contracting states|Kind code of ref document: A1 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
    2020-08-19| AX| Request for extension of the european patent|Extension state: BA ME |
    优先权:
    申请号 | 申请日 | 专利标题
    [返回顶部]